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ABSTRACT 



Instability in liquid monopropellant rocket motors may be 
corrected by the use of a feedback servomechanism. This 
mechanism consists essentially of a pressure pickup which senses 
pressure oscillations in the combustion chamber, an amplifier and 
a variable capacitance in the feed line. It is shown that a feedback 
system with an arbitrary sampling circuit which causes the capaci* 
tance in the line to complete its own cycle of variation once for 
every several cycles of combustion pressure oscillations can be 
made to stabilize the oscillations for all values of combustion time 
lag for a particular motor. It is believed that this system of 
stabilization may be applied to monopropellant motors in general. 
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SYMBOLS AND DEFINITIONS 



time 
time lag 

instantaneous value of pressure 

dimensionless quantity relating pressure and time lag 
steady-state pressure in the combustion chamber 
injector pressure drop 

f/ZAf , injector pressure drop parameter 

distance along feed line between pump outlet and injector 

steady- state mass flow rate 

cross-sectional area- of the feed line 

gas residence time in the combustion chamber 

?rh/zapA©<^ line inertia characteristic parameter 

instantaneous value of pump outlet pressure 

steady-state value of pump outlet pressure 

instantaneous value of pump outlet mass flov/ rate 

- rn I p„-p)/ p^{rrij-w)pump delivery characteristic parameter 

( P"P )/ P » fractional variation of pressure in the combustion 
chamber 

i / 0 <3 , reduced time 

r/Qg , reduced time lag 

/ rVi fractional variation in injector mass flovi^ rate 
reduced amplification coefficient 

reduced angular frequency (real frequency divided by resi- 
dence time) 

K -l- t uo 

instantaneous value of capacitance 
transfer function of sampling circuit 
transfer function of amplifier 



I. INTRODUCTION 
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The lov^ frequency, "chugging", oscillations in liquid rocket 
motors have been attributed to the coupled effects of the feed system 
and combustion chamber dynamics. According to D. F. Gunder and 
D. R. Friant^^\ the pressure oscillations can take place in the 
combustion chamber merely as a result of small changes in the 
velocity of flow in the injected fuel; that is, if for any reason the flow 
rate decreases, the rate of burning and hence the pressure will 
decrease a short time later in the chamber. Then, as the pressure 
drop is felt at the injector the flov; rate will increase and a short time 
later the pressure at the point of combustion will increase resulting 
finally in a decrease in rate of flow from the injection noszle, the 
cycle continuing indefinitely. The explanation for this phenomenon 
involves the concepts of time lag between the instant of injection and 
transformation of the fuel to hot combustion products, and of the 

dependence of the rate of burning on combustion chamber pressure. 

( 2 ) 

M. Summerftsld' ’ has shown how the length of feed line, the 
velocity of the propellant in the feed line, the ratio of feed pressure 
to chamber pressure and the ratio of chamber volume to nozzle area 
affect chugging oscillations and has discussed the trend of changes 
required in the feed system and motor parameters to suppress these 
oscillations. 

( 3 ) 

I/. Crocco' ' has introduced the concept of variable time lag 
basing this concept on a qualitative analysis of the complicated 
processes that affect the propellants from the time of injection to the 
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time of burning. He has formulated a relation between time lag 
and pressure showing that if the pressure is oscillating, the time 
lag is also an oscillating quantity around an average value. Using 
this idea he has demonstrated that self>cxcited oscillations in the 
combustion chamber can exist even in the absence of any variation 
in the injection rate. In the case where the time lag is decreasing 
the burning of the particles that were injected later catch up with the 
burning of those that were injected earlier and an increase in rate of 
burning would result, the opposite being true if the time lag were 
increasing; then, if the variations in time lag coincided with the 
variations in pressure in a certain manner, self>excited oscillations 
would be produced. 

(4) 

H. S. Tsien' ' has demonstrated that a feedback servomech- 
anism may be used to stabilize combustion chamber pressure 
oscillations for all values of time lag, thus making it possible to 
stabilize motors for which a change in configuration parameters alone 
would not correct instability or where such changes would not be 
feasible. The servomechanism Tsien proposed consists essentially 
of a feedback servocontrol which senses the oscillating pressure in 
the combustion chamber and changes the propellant flow by the proper 
amount and in the proper phase to damp the oscillations. Such a 
system is shown schematically on Figure 1. 

A servomechanism in which the variable capacitance is 
required to oscillate at the same frequency as the chamber pressure 
fluctuations is likely to produce at least two important problems 
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arising from the rate at which these oscillations generally occur, 
namely in the order of 100 per second. One is the difficulty of pro' 
ducing sizeable changes in capacitance at such frequencies, and the 
other is cavitation in the fluid. If the variable capacitance consists 
of a piston*cylinder mechanism, it would help to lessen the problems 
to move the piston down on the compression stroke at the same 
frequency as the chamber oscillation and to withdraw the piston 
slowly to the maximum ^capacity position in the time required for 
several pressure oscillations in the chamber to take place, ^’ith the 
proper phasing and amplitude it should be possible to damp out 
pressure oscillations in the chamber. 

The purpose of this investigation is to design a feedback 
circuit that will detect symmetric (say sinusoidal) oscillations in the 
combustion chamber and provide signals to the servomotor that will 
produce the asymmetric motion described above and provide stability 
for all values of combustion time lag. 

As an extension of Tsien's work, F, E, Marble and 

( 5 ) 

D. W. Cox' ’ have shown that a feedback servomechanism may 
stabilize liquid bipropellant motors as well as monopropellant motors. 
Their general method for designing the stabilizing transfer function 
of the feedback loop will be used in this investigation. 



II. DYNAMICS OF THE MO NOPROPEEL. ANT MOTOR 



The basic motor without feedback control considered in this 
investigation essentially consists of a combustion chamber, injector, 
feed line and fuel pump. It tnay be thought of as tv;o dynamic systems; 
one, the combustion chamber and the discharge nozzle: the other, the 
propellant supply system consisting of the propellant pump, propellant 
line and Injector. These two systems are coupled through the require- 
ment that the mass of propellant discharged from the propellant line 
appear in the combustion chamber. The performance of the combustion 
chamber is characterized by a mean pressure p , a mean residence 
time 0g and a mean flow rate no . When the chamber pressure 



varies from p its fractional variation will be denoted (^ = 



P - p 



The performance of the propellant feed system depends upon the cross- 

sectional area A and length P of the feed line, the pressure drop Ep 

across the injector and the local slope o( of the mass flow-pressure 

characteristics of the pump where o< = _ . The 

Po-P<. 

fractional variation in propellant injection rate will be defined as 

m i, — rr> 



The fuel itijected into the chamber does not burn Immediately 
but only after it has remained a certain time t in the chamber. The 
time T required for the transformation of propellant into products 
of combustion may be represented by means of the integral 



4- C p , I q ) = constant 

‘t-T 



where f(p, T ) is a function depending upon the mechanics of the 

O 
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process, and “Tg is the ambient gas temperature in the neighborhood 
of the injector. Since the time lag x depends upon the chamber 
pressure it varies v/ith time in the same manner as the chamber 

(A\ 

pressure does. In fact Tsien'"' has shown that 



where 



dr pet-x) - pCx) 

di p 

^ I 0 C 3 +C p _T^g ) 



r> = 



d lo<^ p 



(3) 

Now, as shown by Crocco' , the continuity of mass flow 
through the combustion chamber may be written as 



d4> 

^ ~r ' • 



(i) 



where S is the reduced time lag 
Furthermore, defining- 



T 






and 




P = 



J 



P 



m 



the dynamic equilibrium of the propellant supply system is given by 




4- 





JJ 4- P<J) ■= O 



( 2 ) 



The presence of ^ and (u In Equations (1) and (2) indicates the 
coupling between the chamber and feed systems. 
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To analyze the stability of the system the following solutions 
may be tried; 

, - 53 C5-i)3 

<pi3) = ^ 4> ^ |a(3)= IVI C3 ) e 



where s= A-»- 'ioo . Using these trial solutions and eliminating 

the common factor Equations (1) and (2) represent a solution 

of the system If the two homogeneous equations below are simul- 
taneously satisfied: 

-si 

( S-t- 1 - n r> c / 5 ® C5 i = o (3) 



P$C5) + J5 -V- CP-v - 4 ) = 



(4) 



For non-trlvlal solutions the determinant of the coefficients of 
Equations (3) and (4) must vanish, that Is 



5-^1 - n -t- n e 



sS 



-si 



- e 



js + iP-^ -^) 



H C'C5>)= O C5; 



This is a transcendental equation in the complex variable s , the 

roots of which determine the stability of pressure oscillations. If 

any root of C'cs)=o possesses a positive real part the system 

is unstable. The Nyquist criterion may be applied to C'C5) to 

(5) 

determine the existence of such roots ' In this method the variable 
5 traces a contour enclosing the right-half complex plane. As s 
moves in the manner shown In Figure 2, that Is, along the imaginary 
axis and a large half circle to the right, the behavior of C'(S) Is 
plotted and the number of complete revolutions noted. This number 
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is the difference between the number of roots and poles of in 

the right half s -plane. If the number of poles can be determined 

independently the number of roots is known and the stability determined. 

-sS 

The term e in C Ci>) makes use of the Nyquist diagram 
difficult. A complete investigation would require plotting the function 
for many values of ^ , the reduced time lag. A different 

technique, proposed by M. Satche^^^ and developed by Tsien^^^, 

requires the function G'cs> to be separated into two parts, the first 

-sS 

consisting only of e and the second of the remaining terms. In 
this case may be expressed as a combination of two deter- 

minants 



£5,(5) = 



5 +- I - n 
P 



Js H- G(P+ h) + , 



I 5) - 



Js - CP+ 



( 6 ) 

( 7 ) 



so that Equation (5) may be written 

-sS ^ 

X), (5) + e X)^( 5 ) ■= o 
(A) - S ^ 

Following Tsien'“' in calling e - S.(s) 
Equation (8) becomes 

9, (^) - 9^(5) = gc5) 



( 8 ) 

£),(s) 

and - — h-< 3 ^cs), 

PjCsi 



(9) 



and ^,(.5) and may be plotted separately. For neutral stability 

Rf(S) -o and becomes a unit circle on the Satche diagram; 

all unstable roots lie within the unit circle. Assuming that GCs ) 
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has no poles in the right-half plane then if either passes 

through the unit circle or encircles the origin the system is 
unstable* >vhile the converse is true. The function G(s). however, 
may have poles in the right*half plane because of possible zeros in 
the denominator determinant . This number of poles may be 

found by applying the Nyquist criterion to and this number 

must be added to whatever number of turns <32. makes about the 
origin. 

As an example of the above criteria applied to a specific case 
a motor with the following parameters may be considered: 



o( = 1.0 J = 1.5 n c 0,6 P=i.O 

Expressing and as follows: 




( 10 ) 




( 11 ) 



the coefficients become: 




( 12 ) 




F or this particular motor 






The characteristic equation for s is 



-sS 



1 + 3 . 1 S -V- 1.5 



o 



2. .5 + o.«^ S 
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For the case where s = l«jo , along the imaginary axis, 

I -V- 3.1 C UO — \.S 
— 

'2-. S“ -V. O.q I Uu 

O-^boj — Z.5 \.3b oui^— 

= _ ^ 

h.zS 0.61 UJ b.ZS -V- o.6l ^ 

The Satche diagram for this function is shown on Figure 3. The 
function g^Ctua) crosses the unit circle and indicates instability 

for values of oo below approximately 0. 9. For large values of 

id 

s off the imaginary axis, let s-Re . Then g 2 .(s) behaves 

l0 y y 

as - Re and as 0 decreases from ir/2 to -rr/Z the 
carve progresses clockwise in a large arc from - oo I to i 
in the left-half plane without encircling the origin indicating no 
further instability. 

References 4 and 5 contain examples of stable and unstable 
motor behavior as determined by analyses similar to the above. 



■!. :| js..,,. 
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m. FEEDBACK STABILIZATION 

The unstable operation of a motor with fixed parameters may be 
corrected by means of a simple feedback servomechanism, consisting 
of a combustion chamber pressure pickup instrument, sampling circuit, 
amplifier, servomotor and a variable capacitance. A feedback system 
in a stability analysis provides a mathematical relation between the 
combustion chamber pressure and the variable capacitance and intro- 
duces a new coupling term in the relations for the simple circuit with- 
out feedback. The relation between the capacitance k ) and the 
pressure may be written symbolically as 



of the amplifier and sampling circuit, respectively. The above rela- 
tion indicates that the line capacitance changes with time in accordance 
with the amplifier signals whose action is controlled by a sampling 
circuit, which interprets the combustion chamber pressure impulses 
in a certain manner. 

A new set of simultaneous differential equations arises from the 

( 5 ) 

introduction of the feedback loop' 

d<b , 



I = Fj. 







where F, — 



and F, — are linear differential-integral operators 




o 
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Using trial solutions 

SE Si Si 

K(i.)=K(s)e ^ cj) (i) = (| ( s) e ^ |uCZ)=M(s)e 

a set of homogeneous simultaneous equations similar to (i) ancj (>4) is 
obtained, 

-sS - si> 

(S + l- n-vne )$(s^-e n/\CS)=-o 

P§CS) + j^Js H- )+ i]ivus^ + -^(P+-^)s] K(S)= O 

F. (S) (S) $ CS) - K (s) = o 



•where F, (S) and F^.(s) are the transfer functions of the sampling 
circuit and the amplifier circuit, respectively. The stability of the 
system now depends on the roots of the third order determinant, 

formed from the coefficients of the above set of equations 

-sS -s S 



S-v i-n-vne -e o 

F,(S) Fits) o -( 



(13) 



This determinant may be conveniently expressed in terms of , 

Equation (6), Equation (7), and a third minor defined as £>4 (s) 



o 



£>4(5) H 



The new characteristic equation obtained by setting the stability 
determinant (13) equal to zero is 



(14) 









-V e 






(s^ k-, (s) F, CS) 



n^(s)] = 



( 15 ) 



For stability anedysis using the Satche diagram Equation (15) may 
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be written 

-s^ D.cm 

e = o ( 16 ) 

Fi(s) Cs) 

where 9 ,(^> and gj,(s) now correspond to the first and second 
terms, respectively, of the above equation and may be plotted 
separately on the Satche diagram. Stable motor operation is to be 
assured by choosing the form of P, cs) and F^cs) so that the 
stability criteria are satisfied. It should be noted that (S) , 
and JD;csi are functions of the fixed motor parameters and are 
independent of the feedback functions. 

Generally, instability in the basic motor arises from the inter- 
section of the g^ci.oo) curve with the unit circle as in Figure 3, 
where the reduced angular frequency is in the order of unity at the 
point of intersection. If for higher frequencies the function 
shows no further instability as is generally the case it is necessary 
only to modify move clear of the unit circle for small 

values of co and to cause the gain of the amplifier to disappear at 
higher frequencies. 
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IV. SAMPLING CIRCUIT FOR FEEDBACK SYSTEM 



The functions of the feedback circuit in this investigation will 
be divided into two parts: one, to provide stabilization and control 
the gain at different frequencies and the other, to transform the 
press ure signal received from the combustion chamber in such a 
manner that the compression stroke of the variable -capacitance will 
be of the same frequency and the return stroke will take considerably 
longer time. A sampling circuit will be responsible for carrying out 
this latter task while the stabilizing function will be carried out by the 
amplifier. 

As a representative example, let the pressure signals received 
by the sampling circuit from the combustion chamber be sinusoidal, 
and the signals which the sampling circuit delivers to the amplifier 
be represented by an asymmetric function, as follows: 

Into the sampling circuit 



sin z 

Out of the sampling circuit 






cos z 



- Cos 






i n4<s rVQl O < ? < TT 

\ r V a I TT i < Zn q 



The input and output functions are shown graphically on Figure 4. 
Over a total interval of time Ztv<^ there are q pressure cycles in 
the combustion chamber for each cycle of motion in the variable 
capacitance. 
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The function 4-(?.) may be expanded in a Fourier series 

CO 

= ) (Q I, cos — \ -»-bsin — 



k = i 



The coefficients are 









Tt 



^ Jo 



ic?.) COS-i^ £ 



h = — 

TTC^ 



f r Lr 

J 4(4) Cos-^i d-i 

o 

Replacing i{z) by the required output functions in the above equations the 

I 

following general expressions for the coefficients are obtained 






1 cos ? cos ' 



p2rrq 



TTc^ 



Jl slo TT 

TT 



J 



COS cos — ^ di 



Z-TT 

Zq-. 



q. 



q k 






- k" 



^ J 



Czs-')^ 

/•■^ f'Zn'q 

cos 2 sio — ^ dt= !— 

o 



Cos — . ^ . 5 i n ? d^; 
Z<i-> 



_k 

TT 



l COS-^ TT 

a 



u" 



k 

cos -^TT 
S 

U.*’ 



I (.'i«i- •')’■ 

The second expressions for and above may be used to 
evaluate the coefficients wheti k .For the ca.se where q= k 
the integral expressions used directly are more convenient. 

Anticipating the use of an illustrative example the coefficients 
for the first eight terms of the Fourier series using a value of <q= 4- 
have been calculated and are tabulated below. 
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k 




b. 


1 


0. 350 


-0.821 


2 


0.226 


-0. 1731 


3 


0. 174 


0.0643 


4 


0. 125 


0 


5 


0.0794 


-0.144 


6 


0. 0420 


-0.0535 


7 


0.0154 


-0.0303 


8 


0 


-0.027 



Now, It is necessary to obtain the transfer function of this 
sampling circuit in terms of the complex variable s eventually to 
investigate the stability of the system. The transfer function is the 
ratio of the I,aplace transforms of the output and input functions of 
the sampling circuit. 

The Laplace transform of the output function is 



,oo 



_ 



-s=e 



S (a^ Cos ^ ^ -t- b sio— ?)e 






-I 



k=. k-. 



K A. ■ I 

k=. ^ 



sV -V 



- (a^qs ^ 



Similarly for the input circuit 



CXS 



- 5 ^ 

slots d€ ^ 






o 
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Therefore, the transfer function of the sampling circuit may 
be written as 

CO ^2. 

^ q( s -.- Ji \ 
ic=i T q»' 

Introducing a term dcs in (17) to take into account the damp- 
ing which exists in a real circuit because of resistance elements the 
function F. (s) will be modified as follows: 

oO , 

^ ^ \ 2cs -+- s 

F,cs) ^modi-f’.ed ° 2 . q ^ ( 18 ) 

It is now possible to proceed to the development of the appro- 
priate amplifier transfer function and to determine the overall 
stability of the system. 
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V. AMPLIFIER CIIARACTERISTICS AND STABILITY ANALYSIS 



An appropriate procedure for determining the amplifier 
characteristics is^®^, (l) to find the approximate behavior of 
for s=itu such that does not intersect the unit circle 

and (2) to determine the function F^c^) so that it has the prescribed 
behavior for small values of co and vanish for large uj . The 
first part of the procedure may be carried out in the following 
manner. Each of the terms in the expression for ^ 

JD, csi 



= - 



JDj (s) -4- P.(5) FjCS) CS) 



(19) 



obtained from Equation (16) may be expressed as a power series in 
Thus 

Pr i (q) I 1 Z. 

<D^(^) - d, + d, s d, s 





w^(s') - d. " d; s 




and 


(S) = B,s -V 0 


(20) 


where 


B,. 5 (p* 4) _ e>,- j 


(21) 



The are obtained by expanding the determinant of Equation (14). 
The terms .0,(5) , jD^Cs) and cs) depend only on the fixed 
parameters of the basic motor. Equation (18) representing the 
transfer function F, cs) of the sampling circuit may be expanded 
as follows 



F 2CS 

2 



s^( 



I 






= ) (i + Zcs+sM 

l< 



\ + 
^7 






b^, ^ 

(zcs -i- 



I - — (zcs-*- S^) -k- — (Zcs -V 



Q^q 

I d- s 

W 
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J 



It is appropriate to choose the following form for the amplifier 
transfer function 






I _L 
s 





(23) 



where the 
now be written 



<3 




are unknown constants. The function 3a 

^ d/^‘ s' ) 

(Ao + A.s* A^5^) -(^’^^6.5+ e,sM 



may 



(24) 



where » for brevity. represent the coefficients of s in the approxi* 
mate expression for F, (s> , Equation (22). The reason for choosing 

the form for F^(s) of Equation (23) becomes clear from inspection 
of Equation (24), since, for s c 0 






d/"' 



d ° ■+ A « 



(-•) 



B. 



and 9a remains finite. Now the requirement that the 

curve lie outside the unit circle for small co may be met as 
follows. Prescribe hy the polynomial 

(i. w) = "So ^ "V a C*- (25) 

choosing the so that stability will be obtained. Then there are 
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two expressions for g^cs) valid for small , namely 

9iCs^ = (26) 

and Equation (24). Since these two expressions are power series 
expansions in s they are equal term by term. The results of setting 
coefficients of like powers of s equal to each other are expressions 
relating the unknown constants and the known coefficients from 

which the values of may be calculated. Thus, after performing 
the necessary operations and retaining terms in s of no higher 
power than s^ the following expressions for the are obtained: 




• 20 - 



The second part of the procedure for determining the complete 
amplifier transfer function Is to obtain an expression which agrees 
with Equation (23) for small values of s and will vanish for large values 
of s so as to leave the original Satche diagram unchanged for large s. 
Consider the following choice of 



Fz(s) = 



1 



(-0 



I d s b s 



I -V h 5 



(30) 



vi^here a, b, and h are constants. For small values of s 

r C-0 

Fz (s) = — — ( t + as + bs^ ) 
s 



and the coefficients a and b must become 



a 






(o> 



(-0 




(31) 



In order that may agree with Equation (23). For large s, 

I 

( 5 ) will behave as 7 , causing the second term In the 

h s 

denominator of In Equation (19)» that Is, the product 

F, ( 5 ) F.jCs^ ^ 4 ( 5 ') , to behave as s® so that the value of 

becomes nearly that of the basic motor without feedback. The value 
of h must be such that stability will be obtained for all values of s 
and will be determined largely by trial and error. 

The exact expression for the function 3 ,^( 5 ) becomes then, 
for all values of s 



g^(s) = 



( o; 



(o) 



d, 5 cl t S ) 



I (o^ iCO 

d. d. s 



oO 

I' 






(^'*'2CS-t-5^)Cqo^.S -4-kb^) (i-^ OS-4- 



Ji- + 2CS-+- 



1 + 



(b.5+ B,s'") 



( 32 ) 
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To illustrate the application of the stabilization procedure 
consider the unstable motor for which the Satche diagram is shown 
on Figure 3. To provide stabilization at low values of to the curve 
of c^^clu~>) must be moved away from unit circle. Choosing the 
following values of 

= -1.25 ^ * *2.0 ^ * “0,25 

and using Equation (25) for the dashed line on Figure 3 is 

obtained. From Equations (27), (28) and (29), using a value of 
c s 0. 1, and summing over the first eight terms of the Fourier 
coefficients, the -(j' are 

0.304 ^ -fj'a 0.032 ^ 5.34 _ 

which fixes the feedback loop for small values of co . To obtain the 
general expression for the amplifier transfer function the constants a 
and b are, from Equations (31), 

a s 3. 395 b a 17.48 

and becomes 

_ 03oA 3.3«HS S 

l-t.O'l = ^ 

s i -V h 

Choosing arbitrarily a value of h = 1 the function becomes 






2.5 -V 



- ( I + 3, I S -V- I. ') 



^-(l-t-O.ZS + S^^4(3tS-*- ) 0.304 U.46s'' 



tl— 



u. 



— O.'Z-S -4- 3 
|6 



\ -h 5 



(i,55 -4- 
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F or the condition u’here s = i uj 



— ( I — 1.5 co^ -V ? . I L OO 1 



c^^(t Uj)- 



^.S-V OS\\oJ 



(i- Uj 4- o.^ I t4j^( W.b^+- uj ) 



U— 



W--I 






— U3 4* Q,Z I <j0 



\.0M 



V4-03 



0-304 - S.iZoj 

ooO + <-«.>'*) 



\(r- l.bUj +l.5ioj) 



This function is plotted on Figure 5 and sho'ws failure to provide 
complete stability since it cuts through the unit circle. 

For a value of h=io however, the function stays 

clear of the unit circle for small values of co , as shown on 
Figure 6 (on this plot also are shown the unstable curves of the motor 
without feedback and the curve for the approximate behavior of 
with feedback for small values of co ). As co increases, 
moves down toward — ool ; as s traces the arc in the right~half 
plane shown on Figure 2, proceeds in a large clockwise arc 

in the left half-plane and returns to the positive imaginary axis and 
finally traces a mirror image of the contour below the real axis. The 
closed contour, therefore, does not enclose the unit circle. 

To check for unconditional stability the Nyquist criterion must 
be applied to the denominator of gz^S) , namely F,C5)fi.(s) , 

As shown on Figure 7, the plot of this function makes no revolutions 
around the origin demonstrating unconditional stability for all values 
of time lag. 

Of special interest is the behavior of in the vicinity of 

CO = 0. 9. For frequencies close to this value the proximity of the 
curve to the unit circle shows that pressure oscillations in the com- 
bustion chamber may be poorly damped. 
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VI. CONCLUDING REMARKS 

The scheme of introducing a sampling circuit in the feedback 
loop to transform the pressure signals in a certain manner and 
alleviate probable mechanical difficulties in the variable capacitance 
does not appear to make the problem of designing a stabilizing 

4 

circuit more difficult than usual. The question now arises as to the 
possibility of designing a practical sampling circuit with real 
components that will transform symmetric signals to non- symmetric 
signals. It is necessary to introduce an actual nonlinear circuit into 
the problem and to analyze its stability and performance. In all 
probability the technique of analysis must be modified. 




I 
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